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ABSTRACT 

Here we present the weak lensing results of A1758, which is known to consist of four subclusters 
undergoing two separate mergers, A1758N and A1758S. Weak lensing results of A1758N agree with 
previous weak lensing results of clusters 1E0657-558 (Bullet cluster) and MACS J0025. 4-1222, whose 
X-ray gas components were found to be largely separated from their clusters' gravitational potentials. 
A1758N has a geometry that is different from previously published mergers in that one of its X-ray 
peaks overlays the corresponding gravitational potential and the other X-ray peak is well separated 
from its cluster's gravitational potential. The weak lensing mass peaks of the two northern clusters 
are separated at the 2.5 cr level. We estimate the combined mass of the clusters in A1758N to be 
2.2 ± 0.5 X IQ^^Mq and r2oo = SSOOt}™ kpc. We also detect seven strong lensing candidates, two of 
which may provide information that would improve the mass measurements of A1758N. 
Subject headings: dark matter — galaxies: clusters: general — galaxies: clusters: individual (Abell 
1758) — gravitational lensing: weak 



1. INTRODUCTION 

Current cosmological models show that baryonic mat- 
ter comprises 4.6% of the total mass distribut ion of the 
universe an d that dark matter comprises 23% ( Komatsu 
et al.|[2010 ), thus the dark matter component is approx- 
imateiy hve times more abundant than baryonic matter. 
Baryonic matter follows the distribution of dark mat- 
ter in the universe and in order to test its existence the 
baryonic matter must be separated from the dark matter 
component, but this is problematic because so far it has 
only been detected during the merging process of galaxy 
clusters. 

Galaxy clusters are permeated by X-ray gas, which 
makes up approxi mately 90% o f the b a ryonic matter in 
the cluster (e.g. David et al. (1990); Gonzalez et al. 

the X-ray gas 



As merging clusters interact 



( [2007| ). 

clouds are slowed through ram pressure while the galaxies 
and the coUisionless dark matter continue nearly uninter- 
rupted (Furlanetto & Loeb 2002 ). Thus, merging clusters 
can provide a unique view where, given the proper con- 
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^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5-26555, 
under program 11194. Also based on data collected at the Sub- 
aru Telescope and partly obtained from SMOKA, which is oper- 
ated by the National Astronomical Observatory of Japan. 

^ Sloan Fellow 



figuration, the dark matter component is separated from 
the X-ray gas. The presence of dark matter can be in- 
ferred by showing an offset of the cluster's overall gravita- 
tional potential from its massive X-ray plasma via gravi- 
tational lensing analysis. Gravitational lensing, which is 
the bending of light of background galaxies as it passes 
near a massive structure, is independent of the type of 
matter present and thus indicates the ove rall gravita- 
tional potent i al of t he cluster (see review by Bartelmann 



(2001))- 



& Schneider 

An ottset between the X-ray gas distribution and the 
mass inferred from gravit ational lensing ha s been de- 
tected in the B ullet cluster dClowe et al. | 2006 [) and MACS 
J0025.4-1222 (Bradac et aL||2008 p. Alternate gravity 
models, which replace the need for dark matter with 
gravi ty being stronger than Newtonian on Mpc scal es 
(e.g. [Milgroml ( p!983| ; |BekenstehI| ( [2004| ) ; |Moffat| ( |2006| ) , 
do not completely model the mass ot the merging clusters 
with the optical and X-ray components alone. Instead, 
alternate gravity models need some smaller amount of 
dark matter compared to cold dark matter models, but 
still require at least twice the amount of matter in 
some dark form as is present in the baryonic component 
( [Brownstein fc Mofrat||2007[ |Feix et al.||2008| ). Ahernate 
gravity models that recover the mass distribution of the 
Bullet cluster can likely model MACS J0025.4-1222 be- 
cause of their similar geometries (the X-ray peaks are 
between the lensing peaks). However, given a system 
with the same mass and different geometry, alternate 
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gravity plus dark matt er models get a diff erent answer 
than general relativity (Angus et al. 20061. Results are 
needed from additional merging clusters with different 
geometries in order to determine which model is correct. 
Additionally, merging clusters are one of the few ways to 
measure the self-interaction cross section of dark matter. 

A1758, which is composed of four clusters, is undergo- 
ing two separate mergers. Evidence that the two clusters 
in A1758N (both are 7 keV clusters) are merging includes 
a projected separation of 800 kpc between its two clus- 
ters, a disturbed X-ray gas distribution, and a separation 
of the X-ray gas from one of the clusters. Evidence that 
the less massive clusters in A1758S (both 5 keV clus- 
ters) are undergoing a merger is that there is very little 
projected separation between the clusters and it has a 
disturbed X-ray morphology (David & Kempner 2004). 
There is no evidence of a merger between A1758iN and 
A1758S in the X-ray sig nature and they have a projected 



Observations of A1758 were carried out on the Subaru 



separation of 2.0 Mpc ( [David fc Kempner 2004 Okabe 
& Umetsu 20081. A1758M introduces a new geometry 



that is different from previously published mergers. One 
weak lensing peak overlaps an X-ray peak while the other 
weak lensin g peak is clearly separat e d from the X-ray 



comp onent (David & Kempner 



2004 Okabe & Umetsu 



20081). The dyna mical state ot A1758K1 given in |David fe 



Kempner (2004) is that the southeast cluster is movmg 
to the sout heast and the north west cluster is moving to 
the north. |David fc Kempn"er] p004 ) conclude that the 
clusters in A1758S appear to be at the point of closest ap- 
proach, but could be a projection effect. However, ^Haines, 
et al. (2009) argue against A1758S being at closest ap- 
proach oased on its lower brightness in the core. There 
is very little redshift info rmation for A1758 (z — 0.279). 
David h K empner (2004) show that A1758N and A1758S 
have a velocity difference of less than 2100 km s^^, which 
is consistent with infall velocity at the redshifts they fit 

from spectral analysis of the X-ray gas^ 

The previous lensing analysis by [Okabe fc Umetsu] 
(2008) was done with two band filters! This analysis 



improves upon their work by using deeper images to get 
better shape information of background galaxies and uses 
a third filter, which allows for better separation of clus- 
ter galaxies from the background. As a result, this work 
has a more detailed weak lensing shear signal and mass 
reconstruction of A1758. 

Throughout this paper, we use the cosmological pa- 
rameters = 70 km sec~^ M pc~^, fijv/ = 0.27 and 
VLty = 0.73 (komatsu et"aLl[2010| . 



2. OBSERVATIONS 

For our analysis, we used observations of A1758 taken 
with two observatories, with the goal of obtaining a broad 
ground-based field of view (Subaru) for a weak lensing 
analysis of the cluster with improved spaced-based imag- 
ing (Hubble) of the cluster cores. Here we describe the 
images we acquired and the reduction process used to 
prepare just the ground-based images for the weak lens- 
ing measurements. In the end, the Hubble images were 
not adequate for getting shape measurements, but were 
used to identify strong lensing candidates. The weak 
lensing analysis is described in Section [3] 

2.1. Subaru 



telescope's Suprime Cam (Miyazaki et al. 2002), which 
is a wide-field camera whose field ot view is 34' x 27'. 
Images in B and V passbands were taken on April 7, 
2008 using time acquired from the Gemini time-share 
program. The weather for the night was photometric. 
Each of our images were 300-sec exposures, which we 
chose in order to minimize the number of saturated stars 
while keeping CCD readout overhead to an acceptable 
level. In total we took 6 images in B and 16 images 
in V, giving 30 and 80 minutes total integration time 
respectively. The telescope was dithered by 20" after 
each image was taken to fill in chip gaps and allow for 
night sky flats to be created. The B images had a range 
of 0'.' 80— 0'.' 91 seeing and a typical value (of the combined 
images) of 0'.'85. The V images ranged from 0'.'76 — 0'.'90 
and had a typical value of 0'.'81. 

To get a third passband, we used archived Subaru 
Suprime Cam Rc images, downloaded from SMOKA 
(Baba et al. 2002) which were taken on July 17, 2007. 
We used 6 images, each with a 500-sec exposure time, 
giving a total exposure time of 50 minutes. The images 
in Rc were taken with a different rotation angle, there- 
fore the field of view of the Rc band did not completely 
overlap that of the B and V images, but did cover both 
the A1758N and A1758S regions. The seeing in the Rc 
images ranged from 0'.'83 — 0'.'92 and had a typical value 
of 0'.'87. After rotating the Rc image into the orientation 
of the B and V images the typical seeing increased to 
0'.'90. 

Image processing was as follows. A master bias file 
was created for each of the ten chips by averaging pixel 
values within 3 ct of the initial mean. These master bias 
files were subtracted from the science frames and then a 
linear fit was made of the overscan strip to subtract any 
residual bias changes. A master flat was flrst created 
for each chip from the bias subtracted dome flat images 
and each chip was smoothed and visually analyzed for 
steep gradients and dust or other undesirable chip fea- 
tures. Typical pixel values were kept in these normalized 
flats in the range of 0.86 — 1.03 and the remaining pixels 
were marked bad. The four corner chips, however, had a 
steep gradient in the far corners of the CCD array and 
pixel values were kept in the range of 0.75 — 1.07 to reject 
the worst parts of the flats while keeping as much of each 
corner chip as possible. The science frames were divided 
by the master flat of the corresponding chip. Another 
set of fiats was created from the fiattened science frames 
and having multiple pointings allowed for excluding in- 
dividual chips with saturated stars. The science frames 
were then divided by these new flats according to chip. 

The sky background was subtracted from each flle by 
making a smoothed image of all local minima within 5 cr 
of the mode of pixel values in order to avoid minima in 
saturated stars. The chip images were scaled so that they 
had equal gain and the same sky levels. 

The brightest stars and galaxies in A1758 were then 
aligned with 633 matching objects in the USNO catalog 
by fitting a P* order polynomial (3 free parameters per 
chip: dx, dy, and rotation) to get the ten chip planes into 
the detector plane. We kept one chip fixed and found a 
solution to the other nine chips in relation to it (27 free 
parameters, which we kept constant for all images). We 
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simultaneously created a 7*'' order 2D polynomial solu- 
tion for each image from the detector plane to the sky 
(36 free parameters per image) with 43,119 total stars 
(stars per image ranged from 521 — 699) from 75 images 
taken on April 7, 2008, which is much larger than the 
(36 X 75) + 27 free parameters needed for the solution. 
We created a separate 7'^ order solution for the Rc im- 
ages taken on July 17, 2007 with an average of 695 stars 
per image. Our mapping to the USNO catalogs was ac- 
curate to 0'.'43 (median), consistent with the accuracy of 
the catalog itself. The consistency of the stars between 
images in this mapping process was accurate to within 
0'.'034 (median) and there were no large scale gradients 
in the residual offsets in these measurements. 

Standard stars taken on April 7, 2008 include the fields 
SA98, SA104, and two pointings of SAllO (the fields 
taken on July 17, 2007 include SA104 and SA107) and 
were reduced in the same manner as above. Standards 
that lay on bad pixels were discarded if the pixels lay 
within several half-light radii of the centroid. The zero- 
point magnitudes were calculated as a function of color 
and airmass with rms variance in _B, V and R of 0.031, 
0.014 and 0.019, respectively, and is the primary error of 
our photometry. There were no large scale gradients in 
the zeropoint levels or errors by position. 

2.2. Hubble Space Telescope 

We obtained observations of the A1758N substructure 
field using a 28 pointing mosaic with the Wide Field 
Planetary Camera 2 (WFPC2) in the F606W filter in 
cycle 16, converting the original 3 passband, 4 pointing 
ACS proposal to a single passband WFPC2 mosaic after 
the ACS shutdown the night after the cycle 16 proposals 
were due. The observations were taken between October 
23rd and 28th, 2007. The mosaic was arranged in two 
layers, a 16 pointing first layer and a 12 pointing second 
layer. The pointings were taken with the same approx- 
imate roll angle, with the second layer pointings offset 
by roughly half a WFPC2 chip in size (~ 40"). The off- 
set was done in order to attempt to identify the charge 
transfer inefficiency (CTI) pattern in the ccds and cor- 
rect it as part of image reduction in a similar manner to 
how the CTI i n the Advance d Camera for Surveys has 
been corrected ( |Massey|2010| . The WFPC2 CTI was so 
bad, however, that we were unable to create a suitable 
correction scheme that would remove the CTI from being 
the dominant systematic error in galaxy shape measure- 
ments. 

To process the images we used the standard calibration 
data and a modified version of custom software from the 
HAGGLes project (Marshall et al. 201 1, in preparation), 
which is based upon MultiDrizzle (Koekemoer et al. 
2002 1 . Specifically, we adapted the code to be compati- 
ble with WFPC2 data and disabled iterative astrometric 
matching due to the limited WFPC2 fleld-of-view. The 
results are the same as using MultiDrizzle. To create a 
final mosaic we utilized the S'C^MP vl.4.6 and SWARP 



lens identification section (Section 4.6) 



v2.17.1 packages (Bertin et al. 2002). The final mosaic 
has a mean internal astrometric uncertainty of 0'.'021. 
Due to the lack of a CTI correction, we found that we 
could not trust the weak lensing measurements enough 
to believe any substructures in the HST data, and so 
have used only the Subaru images for the weak lensing 
analysis below. We use the HST images in the strong 



3. WEAK LENSING ANALYSIS 

After we processed the Subaru im ages and combined 
them by filter, we used SExtractor (Bertin & Arnouts 
1996) on the final images to create a catalog of objects 
and measure their brightness. Then we used a modified 
version of the imcat software]^ to measure size, shape, 
and significance of the objects in the catalog. The resul- 
tant catalogs were purged of objects with less than 10 a 
(to not dilute the shear signal), with an extremely high 
ellipticity (0.6, to reject cosmic rays that are nearly par- 
allel to the detector), or that had any bad pixels (which 
would bias the centroid, fiux and shape measurements). 

We determined the PSF by plotting all objects in ra- 
dius vs. magnitude space and selecting those objects with 
constant radius (~ 2.0 pixels) and V as faint as ^ 24 
mag. These objects are mainly stars, but may include 
some small, faint galaxies. These point-like (or very 
small) objects have shapes that are entirely (or heavily) 
infiuenced by the PSF and provide us a way to model 
it in each image. We want isolated objects so that close 
neighbors do not negatively infiuence our PSF model. 
After removing objects that had neighbors within 30 pix- 
els (~ 15 half-light radii), the remaining objects were 
modeled using a 10*'' order 2D polynomial to iteratively 
fit the objects and reject anomalies. The remaining stars 
were then used to calculate the PSF. 

In order to isolate the background galaxies in the full 
catalog, we removed all objects from the catalog that had 
a smaller half-light radius than the largest stars used 
to create the PSF model (~ 2.2 pixels). Objects were 
also removed if they had an unusually large maximum 
pixel value compared to their overall brightness, which 
indicates something out of the ordinary like a cosmic 
ray or a faulty pixel. The objects that remain in the 
catalog are predominantly the cluster, foreground, and 
background galaxies. 

In order to determine the cluster and foreground galax- 
ies, we plotted simulated galaxies from a catalog provided 
by R. Pello (priv. comm.), which is based on spectral en- 
ergy distribut ion templates of the hyperz photometric 
redshift code ( Bolzonella et al.|2000 ). We plotted galax- 
ies with z < 0.3 in — \/ vs. V — R and compared them 
to the galaxies within 3' of the two BCGs in A1758N. We 
then applied a box cut around the knot of red sequence 
galaxies and the trail of bluer (spiral) galaxies. This box 
cut was applied to our catalog to remove the cluster and 
foreground galaxies, leaving just the background galax- 
ies. 

Next, we discarded any background objects if the sky 
gradient around them was too steep, meaning they were 
too close to a bright star or large galaxy and their shapes 
could be biased. Finally, we discarded galaxies dimmer 
than V ^ 27 and R ^ 26, because they are heavily in- 
fluenced by the PSF and their true shapes are not well 
measured. 

We measured the shear of the resultant catalog of back- 
ground galaxies by calculating the second moment of 
their surface brightness and removed the modeled PSF 



using the KS B method (Kaiser et al. 1995 Luppino & 



Kaiser 1997), which removes the ellipticity of the PSl* 



^ http:/ /www.ifa.hawaii.odu/~kaiser/imcat/ 
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Fig. 1. — The left panel is a gray scale image of the 80-minute V-band image of the A1758N cluster merger (z = 0.279); the two BCGs 
are separated by ~ 3' (800 kpc) on the sky; there could also be a radial component. The blue contours represent the weak lensing mass 
reconstruction made from a background galaxy density of 24.0 galaxies/arcmin^. The outer blue contour begins at k = 0.07 and each 
contour increases in steps of 0.045 up to k = 0.34. The red contours follow the X-ray gas mass based on a 41.5 ksec Chandra exposure. The 
NW cluster's BCG aligns with the X-ray gas and the weak lensing peak. The SE cluster's BCG and weak lensing peak are well separated 
from the X-ray gas, which has a bright peak near the midpoint of the two weak lensing peaks. The two small panels show the smoothed 
red sequence cluster light, which has been weighted by flux. The contours in the upper right panel are the same contours as in the left 
panel. The contours in the lower right panel show the 1 — 3 cr errors on the centroid for the weak lensing analysis. 



from the shear of each object. 

Our weak lensing an alysis of A 1758 was applied to the 
STEP2 data set (Massey et al.||2007j) in order to deter- 
mine the amount of shear recovered in this analysis. The 
data sets psfA and psfC are most like the PSF conditions 
of the Subaru tele scope where the seeing is in the range 
of 0'.'6 - O'.'S. The |Luppino fc Kaiser] ( pQ?] ) method has 
a well known bias that under measures the shear and, 
based on our analysis of the STEP2 data, we determined 
that we recovered 0.893 and 0.875 of the shear in the V 
and Rc images. We divided the shear of each object by 
the respective value in order to correct for this bias. 

To justify our comparison to STEP2, we note that the 
similarities between our images and those simulations in- 
clude pixelation, correlated noise, realistic galaxy shapes 
and a realistic PSF. The only significant differences are 
the potential changes caused by cosmic ray clipping rou- 
tines, the exact form of the correlated noise and PSF 
shapes, and that the PSF was constant in STEP2. In our 
analysis and in STEP2 results, correlated noise does not 
appear to significantly affect shape measurements per- 
formed using KSB-type PSF correction, and any poten- 
tial effect would be significantly smaller than the level of 
random noise in individual cluster mass reconstructions. 
Our co-adding routines were designed to err on the side 
of leaving cosmic rays in the center of bright objects to 
avoid pixel clipping that can alter object shapes, particu- 
larly those of stars, and have been tested on simulations. 
All of these potential effects have been shown to change 
the shears in both STEP2 and similar simulations at the 



few percent level, which is small compared to the random 
errors in the measurements presented in this paper. 

We then placed galaxies into bin regions based on 
size and significance and weighted each object accord- 
ing to 1/rms of the shear of the bin. The shear of ob- 
jects with the lowest SN are thus given the least weight, 
w. We reconstructed the mass of A 1758 by applying a 
Fourier transform to the shear from the V and Rc cat- 
alogs separately and found that the two maps were in 
good agreement with each other and with the BCGs. 
We created our final convergence map by using back- 
ground galaxies that appeared in both the V and Rc cat- 
alogs, calculated the final reduced shear of each galaxy as 
g — (w'^gv + w'^gr) / {Wy + w'^) and the final weights were 
added in quadrature; see the convergence map in Figure 
[l] We computed Tjc„t = 3.21 x 10^ Mf:,/kpc ^ by ta king 
the COSMOS photo-z catalog of [Ilbert et aI] ( |2009D , ap- 
plying the same color and magnitude cuts as we used on 
our background galaxy catalog and taking the harmonic 
mean, which was Zbg = 0.784. 

3.1. X-ray Mass 

We have created an approximate projected gas mass 
map using the Chandra X-ray image in the 0.8-3.5 keV 
energy band extracted from the 41 ks obs e rvatio n (OB- 
SID 2213) analyzed in [ David fc Kempner] ( [2004| ) . For a 
hot cluster such as A1758 [Ig = 9-0 keV, David &: Kemp- 
ner), the X-ray emissivity at photon energies <C Te 
depends very weakly on gas temperature (and its likely 
variations across the cluster). Chandra has a peak of 
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sensitivity at ~ 1 keV, which makes the X-ray surface 
brightness in our energy band a good representation of 
the projected X-ray emission measure, EM cx n^. 

To convert the projected emission measure to the gas 
mass requires knowledge of the three- dimensional cluster 
geometry. Unlike the Bullet cluster (Clowe et al. 2006 



Markevitch et al.||2002[ ), which appears to have a simple 



geometry, the geometry of A1758 is more complicated. In 
the absence of any information or plausible assumptions 
on gas distribution along the line of sight, an approach 
frequently used to get a first-order approximation of a gas 
mass map is simply to take a square root of the X-ray 
brightness. Since on large scales, clusters have spheri- 
cal, centrally concentrated density distributions, here we 
use a slightly more accurate approximation. We first fit- 
ted a spherically-symmetric /3-model to the X-ray radial 
brightness profile and created a projected gas mass map 
that corresponds to that model. This mass map was then 
multiplied by a correction factor {Sx / SpY^"^ , where Sx 
is the real X-ray image and Sp is the /3-model image. 
To normalize the gas mass, we used an Mgas — T rela- 
tion from Vikhlinin et al. (2009; technically, we combined 
their Mtot — T and /gas — T relations, which were derived 
from the C handra gas masse s and X-ray temperatures). 



Though the [Vikhlinin et al.| ( [2009| Mgas - T relation is 



derived for relaxed clusters, hydrodynamic simulations 
indicate that it sh ould not be very different for mergers 
( [Nagai et al.|2007[ ). For the A1758 temperature (T = 9.0 
keV) and z, the relation gives Mgas = l.l±0.28x IO^Mq 
in a sphere of radius rsoo = 1.16 Mpc. We normalized 
the resulting gas mass map to have the projected mass 
in an aperture of the radius 0.65r5oo (which corresponds 
to the Chandra field of view) equal to that of the best- 
fit /3-model with the above mass within the r^QQ sphere. 
The resulting gas mass map is shown in Figure IT] 

Simulations, e.g., by Kravtsov et al. (2006) and Rasia 
et al. (2011), indicate that even the clusters in the middle 
of a major merger, such as A1758, follow the M — T 
relation with a scatter of about 20-25% along the mass 
axis (see, e.g.. Figure 6 in Rasia et al. for quantities 
within r5oo), and the Mgas — Mtot relation is even tighter. 
Other errors in our analysis should be smaller and we 
have assigned a conservative 25% error (68% confidence) 
to the gas masses. 

4. RESULTS OF A1758 

Our WL mass reconstruction of A1758N was created 
by using 20,919 background galaxies that appear in both 
the V and Rc images. The overlapping area between 
the V image and the rotated Rc image is ~ 791 arcmin^ 
and results in a density of 24.0 galaxies/arcmin^ that 
is used for shear measurements. We created the conver- 
gence m ap (see Figure [T]) using the KS93 method (Kaiserl 
fc Squires,1993) , which takes the shear catalog (described 
in Section[5]) and directly inverts it using a Fourier trans- 
form. This produces a 2D projection of the mass along 
our line of sight. We smoothed the mass map enough 
to smooth over small noise features, but not so much 
that our weak lensing peaks would blend together. The 
BCGs of both the SE and NW clusters line up well with 
the weak lensing peaks in our analysis, which are labeled 
Peak 1 and Peak 2, respectively. 

4.1. Cluster Light 



Red sequence cluster galaxies were selected by first 
finding the red sequence of the objects within 3' of the 
two BCGs in a color-magnitude diagram and then se- 
lecting all objects in the catalog with the same color- 
magnitude selection. These cluster galaxies were com- 
pared both by flux and by number density; the results 
were consistent with each other. The side panels in Fig- 
ure [l] show the flux weighted cluster light. The cluster 
light image was also used to select random locations away 
from the cluster galaxies for noise measureme nts o f the 
mass in the bootstrap resamplings (see section 4.2). 



4.2. Errors on the Centroid and Cluster Core Masses 

The errors on the centroids for Peaks 1 and 2 were 
calculated by performing '-^ 3 x 10^ bootstrap resam- 
plings of the background galaxies, recreating the conver- 
gence map, and plotting the nearest centroid (calculated 
with imcat's findpeaks routine) with each resampling (see 
lower right panel of Figure [l]). The results are the green 
and red contours, which show one, two, and three stan- 
dard deviations of the centroids of Peaks 1 and 2. Peak 
1 (2) was found to the NW (SE) of the midpoint 1.19% 
(2.02%) of the time. Figure 1 (lower right plot) shows 
that the separation of Peaks 1 and 2 is significant at the 
level of ^ 2.5 CT. Bootstrapping does not detect system- 
atic errors, such as imparting a false shape due to stellar 
shape. To account for this systematic, we measured the 
error that would occur if the stellar shape were not prop- 
erly removed and thus affect the centroid measurement. 
If left uncorrected, the mass reconstruction based solely 
on the stellar shear field at Peak 1 (2) would shift the 
centroid by a maximum amount of 2" (6"). This is a 
smaller shift than the 68% contour levels in Figure [l] 
The KSB m ethod, however, rem oves the PSF down to a 
few percent ( Massey et al.|2007 1 . Effects for the remain- 
ing unsubtracted portion from the stellar shear field are, 
therefore, negligible. 

The resampled peaks of both Peak 1 and Peak 2 that 
were found past the midpoint (that coincides with an 
X-ray peak) were analyzed to determine why the peak 
finder would completely miss one peak and find it at the 
location of the other peak. The resampled convergence 
maps always showed a significant amount of mass at both 
locations, but in these cases the mass signal steadily in- 
creased from one peak, across the intermediate region, 
and continued increasing toward the other peak. Find- 
ing the centroid away from the expected location was due 
to the linear peak finder and not due to a lack of mass 
in the weak lensing signal at the location of the peak 
in question. The noise in these resamplings occurred in 
such a way as to create a smooth mass increase away 
from the peak being calculated. 

To quantify the amount of mass in these regions we 
averaged the reconstructed mass map pixel values (in 
units of k) in identical, non-overlapping circular aper- 
tures with as large of a radius as possible. The radius 
used for the apertures was 0^663, or 169.7 kpc at the 
redshift of A1758. The mean k within each cutout was 
converted to mass by multiplying by Scrit and the area 
of the aperture. 

When Peaks 1 and 2 were found on the opposite side of 
the midpoint, the average mass of these bootstrap resam- 
plings were low by ~ 2 cr. At the same time, the midpoint 
was ~ 0.5 a higher than its average; in this way the peak 
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TABLE 1 

Mass'' 





Map 


L 


Afap/i 


A/gas 


(Afap - Afgas)/L 




(IQi^ Mq) 


(10" L0) 




(1012 Mq) 


(Az/q/Lq) 


Peak 1 


8.60 ± 1.0 


4.08 


210.7 ± 24.1 


7.39 ± 1.85 


192.6 ± 22.0 


Peak 2 


8.58 ± 1.0 


4.28 


200.6 ± 23.3 


7.64 ± 1.91 


182.7 ± 21.2 


Midpoint 


7.51 ± 1.0 


2.09 


358.5 ± 46.9 


8.33 ± 2.08 


318.7 ± 41.7 


Left 


3.30 ± 1.1 


1.07 


306.8 ± 106.1 


N/A 


N/A 


Right 


4.97 ± 1.1 


1.45 


341.0 ± 72.0 


N/A 


N/A 


Peak 3 


5.72 ± 1.1 


2.74 


208.7 ± 40.4 


N/A 


N/A 


Peak 4 


5.54 ± 1.1 


2.20 


251.5 ± 50.1 


N/A 


N/A 


Random 1 


2.48 ± 1.1 


0.33 








Random 2 


-0.280 ± 1.2 


0.14 








Random 3 


0.197 ± 1.3 


0.14 








Random 4 


0.448 ± 1.3 


0.082 









^ Mass measurements are based on identical, non-overlapping circular apertures centered on the k peaks (and othe r regions) from the weak 
lensing analysis. The k within the aperture radius r = 0'.663 = 169.7 kpc was converted to mass (see Section |4.2| for details). 



finder wandered up the mass structure away from the 
original region. 

The statistics of these points of interest in A1758N, 
as well as four other randomly selected regions (away 
from the red cluster light), show that they are Gaussian 
distributed. The distributions of Peaks 1, 2, and the 
midpoint and that of the Random regions are in good 
agreement with each other. In units of lO^^Af0, the 
mean values of the three peaks are: mp^aki ~ 8.60 ± 1.0, 
mpeak2 = 8.58 ± 1.0, rrimidpoint = 7.51 ± 1.0 (see Table 
fl]). The significance of these masses are 8.7 a, 8.6 a, and 
7.6 a. 

The mean mass at Peak 2 and the midpoint are only 
different by ^ 1 cr, so we tested the mass results for linear 
correlation. Using 3.76 x 10* bootstrap resamplings, we 
found that the mass in Peak 1 (2) was correlated with the 
mass in the midpoint with a coefficient of 0.309 (0.211), 
or 59.9 a (40.9 a) above zero correlation. Peaks 1 and 
2 were correlated with a coefficient of 0.0140, or 2.72 a 
above zero correlation. The mass in Peak 1 (2) was larger 
than that in the midpoint 82.8% (80.8%) of the time and 
there was more mass in Peak 1 compared to Peak 2 50.5% 
of the time. 

These correlations are likely due in part to smooth- 
ing used in the mass map and in part because the large 
aperture cutouts nearly touch; the mass in each resam- 
pling moves with the centroid, pushing mass out of the 
apertures — sometimes toward the midpoint. We also 
compared the mass within smaller apertures (half the 
radius, 84.85 kpc). The linear correlation coefficients be- 
tween Peak 1 (2) and the midpoint were 0.191 (0.130), 
or 37.1 fj (25.3 cr) above zero correlation. Peaks 1 and 
2 were correlated with a coefficient of 0.0188, or 3.6 cr 
above zero correlation. With this smaller aperture, the 
mass in Peak 1 (2) was larger than the mass in the mid- 
point 87.6% (92.9%) of time and Peak 1 was larger than 
Peak 2 33.0% of the time. 

4.3. Cluster Virial Masses 

Now that we have determined that there is a consid- 
erable amount of mass at Peaks 1 and 2, the next mea- 
surement of interest is to determine the masses of the 
clusters in A1758N. Fitting NFW profiles to these clus- 
ters depends on the input parameters for the inner and 
outer radii of the fit. Because the weak lensing approx- 
imation breaks down in the strong lensing region, the 



inner radius should be larger than the strong lensing ra- 
dius. No strong lensing has yet been confirmed in this 
system so we used a value of 100 kpc for the inner radius. 
The outer radius that we used was 3350 kpc, which is the 
distance from the clusters in A1758N to the nearest edge 
of our image. This outer radius is much larger than r2oo 
of both clusters (shown below) . 

We also allowed our color-color selections to vary (see 
Section Is]) because the Subaru filters differ slightly from 
those in ihe R. Pello catalog and because there may be 
systematic errors in our zeropoints. We maximized the 
significance of our cluster fitting by allowing the redder 
side of our box cut to increase redwards by 0.5 mag; other 
shifts in the box edges did not improve the significance 
of the clusters. 

We attempt to simultaneously fit the individual clus- 
ters to obtain their masses, but find that with the current 
data we cannot robustly extract mass measurements. 
Initially, we fit all four clusters simultaneously (two in 
A1758N and two in A1758S) with NFW profiles. How- 
ever, the best fit to one of the weak lensing peaks in 
A1758S was not a well-behaved cluster, but had a very 
large (> 2.0 Mpc) r2oo with extremely low (c < 10^"^) 
concentration. We also tried to fit each cluster in A1758S 
separately (with and without the clusters in A1758N), 
but the fitting behavior in A1758S persisted in these ex- 
treme values for radius and concentration. In order to 
determine the mass in the cluster members of A1758N 
from NFW profiles, we need to decide what to do with 
the members in A1758S. The effect of including both 
southern members is to remove mass from the northern 
members and place it in an extended (and unphysical, 
m > 2 X IO^^Mq) cluster. Excluding both southern 
members greatly increases the radius and mass of the 
clusters in A1758N. The results of fitting just the north- 
ern clusters simultaneously produces different mass re- 
sults than when we fit them one at a time. When fit 
together, one northern member is five times larger than 
the other. When fit separately, the northern clusters are 
approximately the same size. Whether we fit them to- 
gether or separately, the NFW profiles of both of these 
methods have approximately the same significance. Our 
conclusion is that the northern peaks are poorly fit by 
the classic NFW profile and we do not trust this method 
of measuring the cluster component masses. For com- 
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parison, however, fitting two clusters in Ai758N (with a 
single cluster in Ai758S) result in r2oo values (and con- 
centration parameter c) of 1200j^^QQ kpc (c = 13) and 
2000 ± 100 kpc (c = 2.8) with confidence 4.6 ct and 6.8 ct, 
respectively. 

Given the problems with fitting the cores of the clusters 
to measure th eir mass, we also pe rformed aperture mass 
densitometry (Clowe et al. 1998). This routine works 
from the outer radii inward and gives a measurement of 
how much total mass there is along the line of sight out 
to any given radius: 

Cc(ri) = K{r < ri) - R{r2 <r < Vmax) 



(7t) o?lnr 



+ 2(1 - r^Vr^J- W Mdlnr, 

J r2 



(1) 



where R, is the average k within a region and (77) is the 
azimuthally averaged tangential shear. This method as- 
sumes we are measuring the shear 7 and not the reduced 
shear g, but since we are only measuring this at large ra- 
dius, the error from this approximation will be minimal. 
This method, however, gives no indication of the virial 
mass of the cluster, so to determine the virial mass we 



used results by Okabe et al. (2010) from 30 X-ray lumi- 
0.15—0.3 that gives M200 based on 
1.32 (see also 



nous clusters from z 

the mass ratio m2n (< r)/'m^jj{< r) of 
Metzler et al. ( 1999 )). This ratio is due to the additional 
mass along the line of sight and because the cluster mass 
enclosed within a cylindrical volume is larger than the 
3D NFW mass within the same radius. Figure [2] shows 
where this mass ratio equals 1.32 ±lcr, where m2D is 
our aperture densitometry measurement at radius r mul- 
tiplied by Scrit and m^D is the virial mass of an NFW 
profile with r2oo = t- The radii that fall within this range 
are 2170 — 2400 kpc. The ratio m2Dl'mzD — 1-32 occurs 
at a radius of 2300 kpc and using this radius as an esti- 
mate of r2oo we obtain a 3D mass of 2.7 ± 0.4 x IQ^^Mq 
for A1758N. This is higher than th e mass estimate of 
1.6 X lO^^MfD for A1758N given by [David fc Kempner ' 
(2004) (based on X-ray temperature) tor the northern 
clusters, but our measurements using this technique also 
includes a portion the mass of A1758S (2.0 Mpc away) 
within 230 kpc. The mass estimate of 'David fc Kemp 



ner 



(| 2004 p is 1.0 X W-'Mq for the southern cluster. At 
this distance, the amount of A1758S within 2300 kpc of 
the northern cluster is roughly 50% (±10% depending 
on concentration and morphology). This reduces our 3D 
mass estimate of A1758N to 2.2 ± 0.5 x lO^^M©, which 
is still higher than, but in agreement at slight ly higher 
than the l a leve l with, the mass estimate from David fc 
Kempner|(|2004|). 



4.4. Mass-to-Light Ratios 

The M/L ratio in the region of the midpoint in A1758 
is over 300 (see Table fTl), which is high compared to the 
values of our peaks ~zOO. Excess mass in this region 
could be cau sed by a self-interact ion cross section for 
dark matter ( Randall et al. 2008|), so me form of non- 
Newtonian gravity ( |Angus et al.||2007 ), or a higher con- 
centration of galaxies compared to dark matter. To test 
if the high M/L ratio in the area between the clusters 



■ 2D/3D ratio 
Theoretical li 




Fig. 2.— The solid line is the ratio of the 2D/3D mass of A1758N 
and the dashed lines show ztlcr. The dot-dashed, horizontal line 
is the theoretical ratio (~ 1.32) of the 2D/3D masses based on 
measurements of a collectio n of 30 X-ray lumin ous clusters from 
z = 0.15 — 0.3 performed by [Okabe et al. ( |2010^ . The background 
galaxies were put in equally spaced logarittimic bins based on ra- 
dius. Based on this 2D/3D ratio, the r200 of A1758N is 2300^^30 
kpc. 

is significant, we simulated weak lensing observations for 
two NFW clusters the same size as the members of A1758 
and performed the above analysis. We found that the 
smoothing did not preferentially increase the mass in the 
region of the midpoint more than it did in any other di- 
rection. 

We also tested the catalog for the possibility of core 
galaxies being wrongly included with the background 
sample and decrease the reduced shear, and thus the 
mass, in that region. We reduced the number of back- 
ground galaxies in the cores by 10-30% and the overall 
mass was reduced in the cores, but there was no change 
in the ratio of the M/L ratios of the apertures. 

We then compared the M/ L ratios of two other regions 
outside the merger that lie along the line of Peaks 1 and 
2 and are the same distance away as is the midpoint, but 
in opposite directions (see Table [l] "Left" and "Right"). 
These two regions were also founa to have approximately 
the same M/L ratio as the midpoint, which may indicate 
that either there is more dark matter or less red cluster 
light at this distance. 

4.5. A1758S 

We measured the mass of A1758S from the boot- 



strap resarnpling using aperture cutouts (see Section 4.2 
and Table [1]). The masses measured in the two south- 
ern clusters (Peaks 3 and 4, in units of 10^^ Mq) are 
iTLpeaki = 5.72 ± 1.1, and mpeaki = 5.54 ± 1.1, or 5.2 a 
and 5.0 a. 

Peaks 3 and 4 were included in the bootstrap resam- 
pling to determine the errors on the centroid. Of the 
3 X 10^ resamplings. Peaks 3 and 4 were found in the 
region of A1758S approximately 62% of the time. More 
specifically. Peak 3 was found in its own region 18% of the 
time, but 45% of the time it was found in the region of 
Peak 4. The location of the centroid of Peak 4 was found 
in its own region 45% of the time and it was located 16% 
of the time in the region of Peak 3. This is a measure 
of how little mass there is in A1758S and how difficult it 
w as to distinguish between the two cluster members. 

David & Kempner ( 2004 ) find that the southern clus- 
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Fig. 3. — Hu bble Space Telescope ima ge of merging cluster A1758N. Boxes indicate possible strong lensing candidates. Regions A 
match those of jOkabe &: Umetsuj j2b08[ l and F — I are possible new candidates. 



ters are at their c losest approach, but is disputed by 
Haines et al. ( 2009 ) who argue that a lack of star forma- 
tion in the core galaxies of A1758S could be due to them 
being poorer clusters than A1758N. We find that the 
M/L ratio in the southern clusters (see Table fl] Peaks 3 
and 4) match those of the north ern clusters, which does 



not agree with the descri ption of Haines et al. ( 2009 ). 

The mass estimate of Uavid & Kempner ( 2004 ) for 
A1758S, based on X-ray gas mass, is 1.0 x lO^^M©. Com- 
pared to their estimate of A1758N (1.6 x lO^^M©), this 
gives a north/south mass ratio of 1.6. Table [l] shows the 
core masses within identical apertures (r = 169.7 kpc) 
of the northern clusters (Peaks 1 and 2) and southern 
clusters (Peaks 3 and 4). The north/south mass ratio of 
our cluster cutouts is 1.53. 

4.6. Strong lensing candidates 

Our image taken by the Hubble Space Telescope 
(HST) Wide Field Planetary Camera 2 (see Figure p| 
shows possible strong lensing candidates in the core re- 
gion of A1758N. R egions A — E are st rong lensing can- 
didates from ,Okabe fc Umetsu (2008). In our own in- 



dependent analysis we also find candidates D — E (and 
possibly C) and discover regions F ~ I. These are only 
candidates because we do not have spectroscopy or col- 
ors to determine more information about them. Regions 
F — I are not discernible in the Subaru images due to 
ground based seeing. 

In the Subaru Rc images, Regions A and B appeared 
to be candidates, but at the seeing of the Hubble image, 
A is too faint to robustly confirm as an arc and B appears 
to be a faint spi ral galaxy w hose length to width (L/W) 
ratio is ~ 3 (see Meneghetti et al. (2001] for distinguish- 
ing arcs from edge-on spirals). Region U is a borderline 
strong lensing candidate with L/W ~ 5.6, which could 
possibly have another image extending to the west in 
the image, but is too faint to be certain. Regions D and 
E have L/W of ^ 8.4 and ~ 12.5, respectively. Regions 
F—I are newly discovered strong lensing candidates from 
the HST image presented here. Region F is particularly 
interesting as we can see a nearly complete Einstein ring. 
In addition, there seems to be a secondary source be- 
ing lenscd. If such as system were confirmed, and the 
second source is at a different redshift, this would con- 
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strain cosmological parameters like fim and JIa ( Gavazzi 
etan2008). Re gion G is a good candidate, witti L/W 
ot ~ 17.4, tor being a strong arc and for being lensed by 
the edge-on spiral to the east. Region H appears to be 
a strong arc almost certainly influenced by the galaxy, 
but could be a spiral structure. The candidate in Region 
/ has L/W of ~ 6.25 of just the brightest region and 
may extend even further. Most of these candidates ap- 
pear to be primarily caused by individual cluster galaxy 
members, not the cluster core. Only C, and possibly D, 
would provide any information about the cluster masses. 
Without redshifts and multi-color images, we cannot use 
these to improve our mass measurements. With high res- 
olution, multi-color HST images (planned for Cycle 18, 
PI Clowe) , faint sources that cannot be seen in the given 
shallow WFPC2 and ground-based Subaru observations 
will potentially be discovered. 



5. SUMMARY 

A1758 is a region with two merging clusters in its 
northern component, A1758N, and two in its southern 
component, A1758S. The merger A1758N shows a differ- 
ent geometry between its X-ray clouds and weak lensing 
peaks than previously published mergers. Instead of hav- 
ing both X-ray clouds between the weak lensing peaks, 
just one X-ray cloud is between the weak lensing peaks 
and one X-ray cloud overlaps a weak lensing peak. The 
location of the two mass peaks reconstructed from our 
weak lensing analysis are in good agreement with the 
positions of both BCGs of A1758N. 

We calculated the error on the centroids of the two 
mass peaks in A1758N by bootstrapping the background 
galaxy catalog. The two peaks were separated at the 
2.5 tr level compared to there being only one large mass 
with two noise bumps. However, the reason that this 
procedure did not result in a more significant result is 
due to the noise in the mass map reconstructions rather 
than a lack of mass in either peak. We found that 1-2% 
of the time the centroid of one peak wandered across the 
midpoint (or to the location of the other peak) when its 
mass was low by ~ 2 cr and the mass around the midpoint 
was higher than it's average by ~ 0.5 u, creating a ramp 
for the peak finder to walk away from the original peak. 

We performed a multiple cluster fitting of A1758 with 
two clusters in th e north and a single cluster in the south 
(see Section 4.3 for details) and estimate r2oo (and con- 
centration parameter c) for the two northern clusters to 
be 1200i?EJo kpc (c = 13) and 2000 ± 100 kpc (c = 2.8) 



with 4.6 (T and 6.8 a confidence, respectively. However, 



because these parameters are quite degenerate, the fits 
can trade mass between the two peaks and these are 
likely not good estimates of the mass of each cluster. We 
also estimated the virial radius of A1758N by using aper- 
ture mass densitometry, integrating from the outside in. 
This method does not give the virial radius of the cluster 
so we used our 2D mas s at each radius and the ratio given 
by Okabe et al. (2010) that compares m2D/"^3i> = 1-32. 
Our virial radius estimate for A1758N is r2oo = 2300^^30 
kpc and, by the same method, our 3D mass estimate is 
2.2 ± 0.5 X lO^^M©. This is marginally high er than 
the 1.6 X estimate of [David fc Kempner] (2004[ ), 

which they estimated from X-ray temperatures. 

Our analysis does not reveal any new insights about 
the kinematics of clusters in A1758N or A1758S than 



have been previously presented jDavid fc Ke mpner|2004 



Okabe fc Umetsu||2008| [Haines eral.||2009p . The A-ray 
peak that overlaps the N W cluster member may simply 
be a line of sight effect, the merger may be in a later stage 
than the first pass, maybe the NW member is larger than 
the SE cluster and the merger had a large impact param- 
eter, or perhaps the X-ray gas was not stripped from its 
host cluster due to less ram pressure. Getting more red- 
shifts for A1758 could improve our understanding of the 
dynamics of the northern and southern clusters, as would 
improving the S /N of the mass reconstruction using up- 
coming Hubble Space Telescope observations (Cycle 18, 
PI Clowe). 

We detect seven strong lensi ng candidates in A 1 758N , 
three of which coincide with Okabe & Umetsu (20081 
and four new candidates. Wit h the improved se e ing of 
the Hubble, candidate A from Okabe & Umetsu (20081 
appears too faint to discern, candidate B seems to be a 
faint spiral, while C could be a candidate with a multiple 
image, but is too faint to discern. Follow up observations 
with the Hubble Space Telescope are needed to confirm 
the candidacy of these objects and find multiply imaged 
source galaxies in the core regions and help determine 
the strong lensing radius of the clusters. If candidates 
C and D are confirmed, they would help improve mass 
measurements of A1758N. 
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